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Supermassive black hole binaries may exist in the centers of active galactic nu- 
clei like quasars and radio galaxies and mergers between galaxies may result 
in the formation of supermassive binaries during the course of galactic evolu- 
tion. Using the very-long-baseline interferometer, we imaged the radio galaxy 
3C 66B at radio frequencies and found that the unresolved radio core of 3C 
66B shows well-defined elliptical motions with a period of 1.05 ± 0.03 years, 
which provides a direct detection of a supermassive black hole binary. 

The presence of a supermassive black hole binary (or a supermassive binary; SMB) in active 
galactic nuclei (AGNs) has been suggested mainly by periodic optical and radio outbursts (1,2), 
wiggled patterns of compact radio jets (i.e., the indication of precession motions of the radio 
jets) (3,4), and X-shaped morphology of radio lobes (5). AGNs with those characteristics tend 
to be associated with strong radio jets and thus their host galaxies are often giant elliptical 
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galaxies (6). If they were made from mergers between or among nucleated galaxies (7,8), they 
may contain two or more black holes in their central regions because progenitor galaxies may 
contain a supermassive black hole in their nuclei (9,10). Captured black holes will be settled 
under dynamical friction in the core of the merged stellar system and then will form a SMB (11). 
Therefore, it is important to establish the presence of a SMB in an AGN and then investigate 
the true role of a SMB in AGN activity. 

The most direct evidence for a SMB can be obtained by detection of the Kepler orbital 
motion of some emission component close to black holes. One technical problem is that an 
expected separation between the two supermassive black holes (e.g., ~ 10 17 cm) is too small 
to be resolved spatially when we use optical telescope facilities because that separation corre- 
sponds to ~100 micro-arcsecond (//as) if its host galaxy is located at a distance of 100 Mparsec 
(Mpc). However, the use of a technique of phase-referencing very-long-baseline interferometry 
(VLB I) at radio frequencies allows us to achieve such position measurements with the accuracy 
of tens of yuas (12). We looked for the Kepler motion of a radio-emission component in a pair 
of radio sources, 3C 66B and 3C 66A; 3C 66B is a radio galaxy at redshift z = 0.0215 (13), 
whereas 3C 66A is a more distant BL Lac object at z = 0.44 (14). This pairing allows us to 
use this source as the stationary position reference to 3C 66B. They are one of best pairs for 
the phase-referencing VLBI, because they are bright radio sources and their spatial separation 
is only 6 arc min. 

We observed this pair at 2.3 and 8.4 GHz with the very-long-baseline array (VLBA) of the 
National Radio Astronomy Observatory (NRAO) over six epochs between 13 March 2001 and 
14 June 2002. Because the separation between 3C 66A and 3C 66B is smaller than the beam 
widths of the VLBA antennas at 2.3 and 8.4 GHz (26 and 6 arc min, respectively), it is possible 
to observe both sources in each antenna beam simultaneously at these frequencies. 

The radio core position of 3C 66B at each epoch was measured with respect to that of 3C 
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66A. We used the local brightest peak in the map, rather than the Gaussian peak, to determine 
the radio core position. The Gaussian peak may have a larger uncertainty as a result of the 
effects of asymmetrical structure of the sources. 

Time variations of the radio core position in 3C 66B at 2.3 and 8.4 GHz (Fig. 1 and Fig. 
2, respectively) can be fitted by an elliptical motion (Table 1). The reduced chi-square value 
obtained when calculated with the observation error is ~ 0.6 at 2.3 GHz and « 0.05 at 8.4 
GHz (12). Because these values are less than unity, especially at 8.4 GHz, the position error 
might have been overestimated. Assuming the reduced chi-square value of unity, we find that 
the revised position errors, Aa' and AS', are 71 and 60 /ias at 2.3 GHz and 9 and 7 /xas at 8.4 
GHz, respectively (12). Although the orbital periods at 2.3 and 8.4 GHz are nearly the same, 
the major axis at 2.3 GHz is about five times as long as that at 8.4 GHz and the position angle 
of the major axis of the two frequencies differs by ps 24°. 

The averaged period of the fitted core motion is estimated to be 1.05 ± 0.03 years. Although 
this period might be interpreted as due to some geodetic effects related to Earth's orbital motion 
around the sun, such possibilities can be rejected for the following reasons. The annual parallax 
of 3C 66B cannot explain the observed motion, because the distance of 3C 66B (85 Mpc for 
the Hubble constant H = 75 km s _1 Mpc -1 ) leads to the amplitude of the annual parallax of 
only ~ 0.01 /ias. Furthermore, the large difference in the length of the major axis between 2.3 
and 8.4 GHz can not be understood from our knowledge of the annual parallax. Only a possible 
variability of the total electron content (TEC) in the ionosphere may explain this frequency 
dependence, because the excess path due to the ionosphere is proportional to TEC v~ 2 , where 
v is the observing frequency (15). The typical position error due to the ionosphere for the 
separation of 6 arc min is estimated to be 0.005 (q^) /ias, which is negligible (16). Another 
possibility is that the radio core in 3C 66B is subject to microlensing effects by a star in the 
Milky Way galaxy near 3C 66B on the celestial plane (17) and is moved by the annual parallax 
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of the star. We attempted to fit the radio core motion by this microlensing model and obtained 
the reduced chi-square value of 1.6, corresponding to a chi-square probability of ~ 10 %. This 
result indicates that the microlensing model does not explain the sinusoidal core motion as well 
as it explains the elliptical motion of the radio core. 

Because the radio core is located at the root of the jet where the optical depth is unity (18, 
19), the radio core at 2.3 GHz is expected to be located at a greater distance from the central 
engine than that at 8.4 GHz. Thus, a precession motion of the jet provides a natural explanation 
for the observation that the major axis of the radio core orbit is longer at 2.3 GHz than at 8.4 
GHz. Although the orbital major axis should be always perpendicular to the mean jet axis, the 
observed inclination angle between the mean jet axis and the orbital major axis is estimated to 
be 63±9° at 2.3 GHz and 39±6° at 8.4 GHz (calculated on the basis of a 50° position angle 
of the jet, which is estimated from the global jet structure at 2.3 GHz). This fact indicates that 
the observed core motion is not dominated by only the simple precession motion. The most 
plausible explanation for this result is that the elliptical motion of the radio core is attributable 
to the combination of the orbital and precession motions of a SMB in 3C 66B, whereas the 
core motion at 2.3 GHz is dominated by the precession motion of the jet, which at 8.4 GHz is 
dominated by the orbital motion of the SMB (Fig. 3). 

According to Kepler's third law, the observed period can be expressed as, 

P = 2-nG- 1 l 2 r^ 2 {M + mY l '\ (1) 

where where G is the gravitation constant, r is the separation of two orbiting black hole, and M 
and m (M > m) are the masses of the black holes, respectively. Eq. 1 allows us to determine 
the mass density p as a function of only P, 

_ M + m 3tt 
9 = (4/3)7rr3 ~ GP^' ( ) 
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The observed result leads to p « 2 x 10 15 M Q pc -3 (where M Q is the mass of the Sun), which 
is larger than those of massive dark objects in NGC 4258 and Sgr A* (20, 21). Because the 
radio core motion at 8.4 GHz is believed to reflect the orbital motion of the SMB more than 
does that at 2.3 GHz, we regard the length of its major axis as the upper limit of the orbital 
radius of the SMB. Assuming that the less massive black hole emanates the jet, we obtain the 
upper limit of the separation and the mass of the SMB of r max xs 5.4(1 + q) x 10 16 cm and 
M max ~ 4.4(1 + q) 2 x 1O 1O M , respectively, where q = is the mass ratio between the 
objects. We note that these parameters are similar to those of a radio-loud AGN, OJ 287, which 
is believed to be an archetypical AGN with a SMB (22). This conclusion is consistent with 3C 
66B's being a giant elliptical galaxy which suggests a galaxy merger. 

The SMB in 3C 66B will merge into one as a result of the gravitational radiation loss. The 
lifetime in this evolution phase can be estimated by the relation in Eq. 3 (23). 



If we adopt q = 0.1, we obtain r max = 5.9 x 10 16 cm and M max = 5.4 x 1O 1O M . These 
values give us t GR ~ 5 yr for the SMB in 3C 66B. This time scale seems to be extremely short 
compared with a whole life time of a SMB system (13). A more reasonable lifetime would be 
obtained for smaller values of q and/or r. Combining Eq. 1 and Eq. 3, we find t GR oc r~ 5 . If the 
true r is a quarter of the corresponding value for the upper limit, a lifetime ~ 10 3 times as long 
is obtained. The expected maximum amplitude of gravitational radiation from 3C 66B will be 
(h) ~ 1CT 11 just before the two black holes merge into one (24, 25). 
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Table 1: Parameters of the fitted orbital motion of the radio core in 3C 66B. The errors in the 
parameters correspond to a change of 1 in chi- square from the value at the best fitted parameter. 
The offset angle is the angle between the major axis and the point where the orbital motion 
started. R.A. center is the relative position in right ascension. Dec. center is the relative position 
in declination. 







2.3 GHz 


8.4 GHz 


Major axis 


(//as) 


243±30 


45±4 


Axial ratio 




0.31±0.17 


0.24±0.14 


Orbital period 


(yr) 


1.10±0.06 


1.02±0.04 


Position angle 


(°) 


113±9 


89±6 


Offset Angle 


(°) 


60±7 


101±5 


R.A. center 


(mas) 


1.441 ±0.020 


0.970±0.002 


Dec. center 


(mas) 


-0.888±0.033 


-1.861 ±0.004 



Fig. 1. Orbital fit calculations applied to the data of the position change of the core at 2.3 GHz. 
Observations were carried out on 13 March, 25 June, and 9 November 2001, and 8 February, 
21 February, and 14 June 2002. We show (A) spatial distribution, (B) time evolution toward 
right ascension direction, and (C) time evolution toward declination direction. The error bars 
indicate the revised position errors suggested from the chi-square value; 71 and 60 //as at 2.3 
GHz and 9 and 7 //as at 8.4 GHz, in right ascension and declination, respectively (12). 
Fig. 2 Same as Fig. 1, but at 8.4 GHz. 

Fig. 3 A schematic view of the geometrical relation between the radio core motion and the SMB. 
Because the orbital velocity of the SMB is added to the intrinsic jet velocity, the jet precession 
can be induced with the same period as that of the SMB orbital motion (3). Therefore, if 
the radio core is far enough from the SMB, the orbital major axis of the radio core should be 
observed as perpendicular to the mean jet axis, as is the normal precession. However, if the 
radio core is very close to the SMB, it should be observed as almost parallel to the orbital 
plane of the SMB. Therefore, the position angle of the radio core orbit depends on the distance 
between the radio core and the SMB, if the jet direction misaligns with the angular momentum 
vector of the SMB orbit. 
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Materials and Methods 

Phase-referencing VLBI 

In conventional observations with VLBI, the correlated phase is affected by 
fluctuations due to water vapor in the troposphere and the fluctuations strongly 
limit the quality and the dynamic range of images. The most effective method 
to eliminate the tropospheric fluctuations can be made by observing simul- 
taneously a target source and another radio source (e.g., a different quasar) 
which is close to the target on the celestial plane. This procedure is called as 
the phase- referencing technique in VLBI observations (S1,S2). In principal, 
assuming that the atmospheric phase fluctuation is almost the same between 
the two sources, we eliminate it by subtracting the phase of the target from 
that of the near radio source. 

In our observations, the data of 3C 66A were first reduced by using the 
standard hybrid mapping method and self-calibration procedures with the 
minimum integration time of 6 seconds available on the NRAO AIPS. We 
thus obtained antenna-based solutions including the atmospheric phase fluc- 
tuation. Then these solutions were directly applied to the observations of 3C 
66B and the final phase- referencing map of 3C 66B was obtained (Fig. SI). 

Error analysis 

The dominant uncertainty of the position measurement in the phase-referencing 
observation comes from unmodeled source structure in the map. In order to 
examine this uncertainty, we compare the local brightest peak with the cen- 
troid of the strong CLEAN components (more than 25 % of the maximum one) 
(S3,SJ f ). This procedure gives nominally us the observational position errors, 
Act in right ascension and AS in declination, 91 fias and 76 fias at 2.3 GHz 
and 40 fias and 28 fias at 8.4 GHz, respectively. 

On the other hand, we derived the standard phase errors from theoretical 
estimates based on models of the propagation medium and the geometry of 
the array, e.g., tropospheric and ionospheric effects, uncertainties of antenna 
positions, the time difference between atomic time and universal time, and 
the polar motion or the nutation of the Earth. The misidentification of the 
reference point in maps due to finite signal-to-noise ratio also contributes to 
the position error. Assuming the typical weather and geodetic conditions in 
VLBA observations, we calculate the corresponding position error in the si- 
multaneous observations of the source pair with the separation of 6 arcminute. 
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The standard position errors, Act in right ascension and AS in declination, are 
estimated to be 10 fias and 16 fias at 2.3 GHz and 7 fias and 12 fias at 8.4 
GHz, respectively. 

As shown in the text, assuming that the reduced % 2 value for the core 
motion fitting is unity, we find that the revised position errors, Ac/ and AS', 
are 71 //as and 60 //as at 2.3 GHz and 9 //as and 7 //as at 8.4 GHz, respectively. 
The revised errors at 8.4 GHz is comparable to the theoretical standard errors, 
suggesting that the uncertainty due to complex source structure is excluded in 
the data at 8.4 GHz. In contrast, the position measurement at 2.3 GHz seems 
to contain effect of unresolved jet structure in the radio core component. 
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Fig. SI Phase-referencing images of 3C 66A at (A) 2.3 GHz and (B) 8.4 
GHz and 3C 66B at (C) 2.3 GHz and (D) 8.4 GHz, obtained at the 5th epoch 
(21 February 2001). The contours are (1, 2, 4, 8, 16, 32, 64, 128, 256, 512)x3cr 
errors. The synthesized beam is plotted at the bottom left of each image, 
with an approximately 5.4x3.3 mas, P. A. =—8° at 2.3 GHz and 1.3x0.8 mas, 
P.A. = -8° at 8.4 GHz. 



2 



References 



51. W. Alef, in NATO ASI Ser., Very Long Baseline Interferometry, Tech- 
niques and Applications, eds. M. Felli & R. E. Spencer (The Nether- 
lands: Kluwer Academic Publishers, 1989), p 261 

52. A. J. Beasley, & J. E. Conway, in ASP Conf. Ser. 82, Very Long Base- 
line Interferometry and the VLBA, eds. J. A. Zensus, P. J. Diamond, 
& P. J. Napier (San Francisco, 1995), p 328 

53. J. C. Guirado, et al., Astron. Astrophys. 293, 613 (1995). 

54. M. J. Rioja, R. W. Porcas, Astron. Astrophys. 355, 552 (2000). 



3 



